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INTRODUCTION 


Quite  often  the  aerodynamicist,  when  he  works  with  or  near  a design  group,  is 
faced  with  the  task  of  estimating  such  important  design  parameters  as  range,  static 
margin,  maneuverability,  dynamic  stability,  time  to  hall  amplitude,  etc.,  tor  a given 
configuration.  Once  the  design  engineer  obtains  the  data  he  is  seeking,  the  next 
question  the  aerodynamicist  must  answer  is  “how  can  1 improve  the  configuration  so 
its  aerodynamic  properties  are  better”?  Ultimately,  an  iteration  cycle  will  probably 
be  made  in  which  several  different  configurations  will  be  considered  before  the  two 
or  three  most  optimum  candidates  are  chosen  for  further  study.  The  important 
point  to  be  made  here  is  that  for  each  of  the  above  configurations,  static 
aerodynamics  (lift,  drag,  and  pitching  moment),  and  dynamic  aerodynamics  (roll 
damping,  pitch  damping,  and  magnus  moment)  must  be  estimated  before  questions 
concerning  such  things  as  range,  maneuverability  and  flying  qualities  can  be 
addressed. 

To  obtain  the  above  set  of  aerodynamic  coefficients,  the  engineer  can  go  one 
of  three  directions:  he  can  conduct  wind  tunnel  tests  which  will  be  costly  and  time 
consuming  and  probably  produce  results  which  are  more  accurate  than  warranted  for 
preliminary  and  intermediate  designs,  he  can  perform  hand  calculations  using 
handbook  techniques'  and  applicable  experimental  data  but  not  have  a good 
accuracy  estimate  of  the  results,  or  he  can  develop  a computer  program  based  on 
analytical  techniques  which  is  efficient  and  produces  accuracies  on  the  order  of 
±10  percent  for  most  configurations.  The  latter  alternative,  although  being  more 
costly  and  time  consuming  initially,  is  the  best  approach  for  long  term  use  and  is 
the  procedure  which  is  followed  in  this  work. 

To  be  of  practical  use  for  tactical  weapons,  the  theory  must  compute 
aerodynamics  for  the  Mach  number  and  angle  of  attack  range  of  most  missiles  and 
projectiles,  that  is,  0 < <3  and  0 < a < 15°,  respectively.  Also,  quite  general 

body  and  wing  geometries  must  be  considered.  This  arises  from  the  fact  that 
projectile  noses  may  be  pointed,  truncated,  or  spherically  blunt.  Another  contributing 
factor  to  the  complex  projectile  geometries  is  the  high  setback  forces  at  launch 
which  means  the  wings  and  canards  must  be  quite  thick  to  survive  the  initial  high 
“g”  loads.  Moreover,  there  may  be  two  ogives  on  the  nose  (in  the  case  of  a fuze 
or  blunt  seeker)  and  a boattail  present  for  drag  reduction  purposes.  Figure  1 
illustrates  the  general  type  of  tactical  weapon  geometries  that  are  encountered  on 
nonairbreathing  configurations.  By  designing  a computer  program  to  handle  such 
complex  geometries  most  missile  and  projectile  configurations  can  be  considered. 
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I.  Guided  Missile/Projectile  (IMonairbreathing) 


Geometry 


a)  Body  — Nose  can  be  pointed,  blunt,  or  truncated  with  discontinuities  along 
surface,  with  boattail  or  flare  allowed. 


b)  Wing  — Leading  and  trailing  edge  can  be  sharp  or  blunt  with  general 
modified  double  wedge  or  biconvex  airfoil  sections  and  up  to  two  sets  of 
lifting  surfaces  allowed. 


II.  Guided  Missile  (Airbreathing) 


Geometry  same  as  for  nonairbreathing  case  except  inlets  added  to  body  or 

wing. 


III.  Spin  Stabilized  Projectiles,  Bombs,  Rockets,  Other  Ordnance 

Geometry  no  more  general  than  that  already  considered  for  guided  missiles. 


Figure  1.  Tactical  Weapon  Configuration  Requirements 


Several  works  existed  previously  which  could  fulfill  portions  of  the  present  goal 
of  a general  aerodynamic  prediction  program,  but  none  of  which  was  satisfactory  in 
its  entirety.  The  most  notable  of  these  is  that  due  to  Woodward.’  Woodward  uses 
perturbation  theory  to  compute  the  pressure  distribution  on  wing-bodv  combinations 
in  subsonic  and  supersonic  flow.  However,  the  bodies  must  be  pointed  and  the  wing 
leading  edge  sharp.  Also,  he  does  not  calculate  the  base  and  skin-friction  drag  or 
the  nonlinear  angle  of  attack  effects.  Moreover,  no  consideration  is  given  to 
transonic  flow  or  to  the  computation  of  the  dynamic  derivatives.  Consideration  was 
given  to  an  approach  such  as  that  followed  by  Woodward  and  others  for  use  in  the 
present  analysis  to  calculate  static  aerodynamics  on  tactical  weapons.  However,  the 


panel  methods,  which  are  currently  in  use  by  many  scientists  in  the  field,  tend  to 
work  quite  well  for  predicting  lift  and  pitching  moment  but  have  limited  accuracy 
in  predicting  drag  for  Mach  numbers  up  to  2.5  or  3.0  and  cone  angles  to  25°. 
Since  drag  is  at  least  as  important  as  lift  for  most  tactical  weapon  performance 
characteristics,  a different  approach  was  followed  for  wing-body  static  aerodynamics. 

Another  method  available  for  calculating  aerodynamics  on  wing-body 

configurations  is  that  of  Saffell,  et  al.3  This  procedure  computes  static  aerodynamics 
on  low  aspect  ratio  missile  configurations.  Its  applicability  to  general  aspect  ratio 
configurations  is  thus  questionable,  particularly  at  small  a.  Furthermore,  drag  was 
calculated  by  handbook  techniques1  and  is  also  quite  inaccurate  at  small  a.  Again, 
no  attention  was  given  to  the  dynamic  derivatives. 

The  other  method  of  practical  use  in  projectile  work  is  the  all  empirical  (IF 
"Spinner”  program.4  This  program  computes  all  six  aerodynamic  coefficients  so  a 
dynamic  stability  analysis  can  be  conducted  on  a spin-stabilized  projectile  without 
going  to  the  wind  tunnel.  Its  short  comings  are  its  lack  of  attention  to  guided 
weapons  and  limited  use  in  body  alone  design  due  to  its  empirical  nature.  That  is. 
effects  of  nose  bluntness,  ogive  shape,  boattail  shape,  or  other  discontinuites  in  slope 
cannot  be  estimated  in  detail.  On  the  other  hand,  it  does  give  total  force 
predictions  which  are  reasonably  accurate  on  spin  stabilized  rounds  which  are  similar 
in  shape  to  the  empirical  data  base. 

To  summarize  the  above  state-of-the-art  in  aerodynamic  prediction,  it  is  fair  to 
say  there  are  several  analytical  or  empirical  methods  available  to  compute 
particular  force  or  moment  in  a given  Mach  number  regime.  However,  these  methods 
have  not  been  combined  together  to  form  a useful  design  tool,  allowing  one  to  not 
only  estimate  aerodynamics  accurately,  but  to  conduct  static  and  dynamic  stability 
analysis  or  trajectory  simulations  on  a given  configuration  over  the  Mach  number 
range  from  I to  3 and  angle  of  attack  range  from  0 to  15°.  It  is  the  purpose  of 
the  present  work  to  provide  such  a tool. 


ANALYSIS 


STATIC  AERODYNAMICS 

The  static  aerodynamics  have  been  computed  previously  and  the  reader  is 
referred  to  References  5 through  9 for  the  details  of  that  portion  ol  the  present 
work.  Only  a synopsis  of  that  work  will  be  included  here  for  continuity  purposes 


Body  Alone  Aeroth  mimics 

A summary  of  the  various  methods  tor  computing  body  alone  aerodynamics 
appears  in  Figure  2.  All  the  methods  are  standard  in  the  literature  (References  10 
through  15)  with  the  exception  of  the  empirical  schemes  derived  for  transonic  lit t 
and  wave  drag  and  the  combined  Newtonian-perturbation  theory  for  calculating  nose 
wave  drag  in  supersonic  flow.  The  combined  Newtonian-perturbation  theory  was 
developed  so  reasonable  results  for  static  aerodynamics  could  be  obtained  at  low 
supersonic  Mach  numbers  for  blunt  nosed  configurations.  Comparisons  with 
experimental  data  indicate  this  method  accurately  predicts  pressure  coefficients  and 
total  force  coefficients  down  to  supersonic  Mach  numbers  of  1.2. 


\ MACH 
\ NUMBER 
COMPONENT  \ REGION 

SUBSONIC 

TRANSONIC 

SUPERSONIC 

NOSE  WAVE  DRAG 

Wu  and  A0Y0MA 
PLUS 

EMPIRICAL 

2nc*  ORDER  VAN 
DYKE  PLUS 
MODIFIED 
NEWTONIAN 

BOAT  TAIL  WAVE  DRAG 

Wu  and  AOYOMA 

2nd  ORDER  VAN 
DYKE 

SKIN  FRICTION  DRAG 

VAN  DRIEST  II 

BASE  DRAG 

EMPIRICAL 

INVISCID  LIFT  and 
PITCHING  MOMENT 

EMPIRICAL 

Wu  and  AOYOMA 
PLUS  EMPIRICAL 

TSIEN  Ist ORDER 
CROSSFLOW 

VISCOUS  LIFT  and 
PITCHING  MOMENT 

ALLEN  and  PERKINS  CROSSFLOW 

Figure  2.  Methods  Used  to  Compute  Body  Alone  Aerodynamics 
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W inn  and  I n terfercncc  Acrody mimics 

The  methods  used  to  compute  the  wing  alone  and  interference  aerodynamics 
are  listed  in  Figure  3.  Again  the  methods  are  standard  in  the  literature 
(References  16  through  22)  except  for  the  empirical  techniques  used  for  wing-body 
interference,  trailing  edge  separation  drag,  and  body  base  pressure  drag  caused  by 
tail  fins  and  except  for  the  theoretical  computation  of  wing  wave  drag  in  supersonic 
flow.  For  the  details  of  these  techniques,  the  reader  is  again  referred  to  either 
References  7 or  9. 


\mach  number 

COMPONENT^.  REGION 

SUBSONIC 

TRANSONIC 

SUPERSONIC 

INVISCID  LIFT  AND 
PITCHING  MOMENT 

LIFTING  SURFACE 
THEORY 

EMPIRICAL 

LINEAR  THEORY 

WING- BODY  INTERFERENCE 

SLENDER  BODY  THEORY  AND 
EMPIRICAL 

LINEAR  THEORY. 
SLENDER  BODY 
THEORY  8 EMPIRICAL 

WING-TAIL  INTERFERENCE 

LINE  VORTEX  THEORY 

WAVE  DRAG 

EMPIRICAL 

LINEAR  THEORY  + 
MODIFIED  NEWTONIAN 

SKIN  FRICTION  DRAG 

VAN  DRIEST 

TRAILING  EDGE  SEPARATION 
DRAG 

EMPIRICAL 

BODY  BASE  PRESSURE  DRAG 
CAUSED  BY  TAIL  FINS 

EMPIRICAL 

Figure  3.  Methods  Used  to  Compute  Static  Wing  Alone  and 
Interference  Aerodynamics 
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ROLL  DAMPING  MOMLN1 

The  body  alone  roll  damping  moment  is  estimated  using  the  same  procedure  as 
tlie  “Spinner”  program.4  That  is.  if  one  knows *C\  as  a function  ol  Mach  number 
for  a given  configuration,  then 


where  the  subscript  represents  known  data.  This  estimate  should  prove  reasonable  so 
long  as  the  data  available  is  for  a configuration  of  similar  dimensions  to  the 
configuration  for  which  data  is  required. 

I or  a wing-body  configuration,  the  wing  completely  dominates  the  roll  damping 
moment  so  it  is  of  primary  concern  to  have  a good  estimate  of  the  wing  alone  roll 
damping. 

For  supersonic  and  subsonic  How  calculations,  the  wing  is  assumed  to  be  thin 
so  that  small  perturbation  theory  can  be  employed  for  How-field  calculations.  In 
addition,  the  wing  is  assumed  to  have  zero  camber  with  aeroelastie  effects  being 
small.  In  transonic  tlow,  empirical  methods  are  used  for  roll  damping  calculations. 
Body  interference  is  estimated  using  slender  body  theory.  Tail  fin  to  tail  fin 
interference  is  zero  unless  the  non-dimensional  spin  pd/2Vm  is  high  enough  so  the 
shed  vortices  impinge  on  an  adjacent  fin. 

The  individual  methods  used  for  calculating  the  wing  roll  damping  moment 
coefficient  derivatives  in  subsonic,  supersonic,  and  transonic  flow,  along  with  the 
interference  effects,  are  discussed  below. 

Subsonic  Flow  (M^  < Mcrjt) 


The  small  perturbation  equation  for  three-dimensional  steady  tlow  is:lh 


*0*0  ^0^0 


O'O 


= 0 


(I) 
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where  the  subscripts  x(),  y0.  z0  indicate  partial  differentiation  and  where  the 
velocity  potential,  0.  is  related  to  the  perturbation  velocities  by 


0,  = u/V  , </>,.  = v/V  , 0,  = w/V 

**0  - Vo  “ *n 


(2» 


The  boundary  conditions  are  that  the  flow  is  undisturbed  at  an  infinite  distance 
from  the  surface,  is  tangential  to  the  surface,  and  is  continuous  at  the  wing  trailing 
edge.  Mathematically,  these  boundary  conditions  can  be  treated  in  respective  order 

by: 


u 


v = w = 0 as  J x2  + y2  + /2  ~+  °° 

c J c c 


(3a) 


vm  • 

m 

VI'  = 0 

(3b) 

uu^xoVi;’yo,zo  1 

= UC  l(X0  )TE  'y0’Z0  1 

vul(xo  *t  i:  ’Vo^o1 

= v8[(x0)T|,y0,z0l 

(3c) 

wJ(xo  *n  -y0’zoi 

= we[(x0)T|:,y0,z0l 

The  boundary  condition  (3c)  is  simply  a definition  of  the  Kutta  condition  which 
requires  the  velocity  at  the  trailing  edge  to  be  continuous  from  the  upper  to  lower 
surfaces.  If  the  equation  of  the  surface  is  defined  as  zQ  = F(xQ,y0),  then  (3b) 
implies 


w<  x()  ,y0 ) 


* K 


9Hxu,y0) 

3xn 


a(x0,yc 


(4) 
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The  angle  of  attack,  cv(  x(|  ,y(| ) can  be  expressed  as 


a(x0,y0)  = 


«n  + py(1/v. 


+ (l(x0- 


x ,.)/V 

ref " • 


(5) 


where  a()  is  the  angle  of  attack  of  the  wing  planfonn,  py0/V^  is  the  induced  angle 
of  attack  at  a point  y()  on  the  wing  span  due  to  a steady  roll  rate  p.  and 
q(  \()  - xrL.f)/V,„  's  the  induced  angle  of  attack  at  a point  x()  along  the  wing  chord 
due  to  a constant  pitch  rate  q. 

liquation  (1)  is  a linear  partial  differential  equation  so  that  superposition  of 
solutions  is  valid.  Therefore,  for  roll  damping  calculations,  one  can  set  a(|  = q = 0 
in  liquation  (5)  so  the  boundary  condition  of  interest  is 


0,o  * py0/v 


(6) 


s’  Equation  (1)  can  be  simplified  somewhat  by  using  C.othert’s  extension  to  the 
Prandtl  Cdauert  transformation.  This  transformation  is  equivalent  to 


x = xo/0,  i 

y = y0  f 

z = Z0  1 

0( x.y ,/.)  = 0(x(J,yn,zo)  ' (7) 


where 


Pi 


= -JT 


M- 


liquation  ( 1 ) then  simplifies  to 


<j)  + <p  + <p  - 0 (8) 

l b us,  the  compressible  How  on  a wing  of  given  geometry  can  be  solved  for  by 
affinely  relating  the  wing  to  another  wing  with  the  properties  of  Equation  (7),  and 
solving  for  the  How  field  on  the  new  wing  at  = 0.  Once  this  is  done,  the 
pressure  coefficient  at  any  point  on  the  wing  is: 1,1 
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('p(Xo-y0>/r»)  = 


(9) 


The  solution  to  I tj nation  (8)  can  be  shown  to  be:17 

d x , d y , 1 10) 

where  AC„  = CP  - C\>  . It  is  required  to  determine  the  pressure  loading  AC,,  over 
the  entire  surface,  hollowing  Chadwick,  Kquation  (10)  is  lirst  ditlerentiated  with 
respect  to  z and  the  limit  as  y.-*0  taken.  The  result  is  then  equated  to  the 
boundary  condition,  Kquation  (6).  to  obtain 


, , x xi 

JJ  (y  - y.  )2  + z2  X 

1 ■>  / 1 
*J(x  - x.  + (y  - y.  r + r 

py 

v^ 


i 

Hit 


AC  p ( x | ,y  j ) 
(y  - y,  )2 


J (x 


x -x, 


x,  r + (y 


y,  r 


d x | dy | 


(II) 


The  cross  on  the  y,  integral  indicates  a singularity  at  y = y, , in  which  case 
Manglers  principal-value  technique16  can  be  applied.  The  details  ot  the  solution  ol 
the  integral  Kquation  (10)  for  ACp(x,y)  will  not  be  repeated  here  as  they  are  given 
in  many  references  (see  for  example.  Reference  17).  Worthy  of  note,  however,  is 
the  fact  that  Kquation  (11)  is  an  integral  equation  for  which  the  wing  loading  AC,, 
is  to  be  found  as  a linear  function  of  the  induced  angle  ot  attack.  py/V^. 


Once  the  span  loading  ACp(x.y),  due  to  roll,  is  known  over  the  entire  wing, 
the  local  rolling  moment  coefficient  at  a given  spanwise  airfoil  section,  y.  is 

x T 1 

(ACp)pdx  (12) 

LI 

where  the  subscript  P indicates  the  loading  due  to  a roll  rate  p.  The  total  rolling 
moment  on  the  entire  wing  is  then 


2 r"'2 

, = — j cc.dy 

ret 


(13) 
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Assuming  the  rolling  moment  is  a linear  function  of  roll  rate,  the  roll  damping 
moment  coefficient  derivative  is 


( , 


(14) 


Supersonic  Flow  (M  > 1.2) 


For  supersonic  flow  past  thin  wings.  I quation  ( I ) is  still  applicable  along  with 
the  associated  boundary  condition  (3b).  Since  the  flow  is  supersonic,  disturbances  in 
the  flow-field  are  not  felt  upstream  of  the  point  of  disturbance.  Ilius,  the  boundary 
condition  (3a)  can  be  modified  to. 


u(0  ,y./.)  = v(0  ,y./)  = w(0  .y./)  = 0 


(15) 


where  it  is  assumed  the  disturbance  occurs  at  x = 0.  It  will  be  assumed  the  wing 
trailing  edge  is  supersonic  (Mach  number  normal  to  wing  trailing  edge  is  greater  than 
one)  so  that  the  Kutta  condition  need  not  be  applied.  The  solution  to  liquation  ( 1 1 


'o  o 0 1 


yo-yo)  . . 

~ r"-'  ' j ^ <,  dx.dy. 

J(x(  - x )‘  - (]2(yQ  - y,  )3 


(16) 


where  S(  is  the  wing  planform.  The  lifting  pressure  coefficient,  due  to  a constant 
roll  rate  p,  is  then  related  to  the  velocity  potential  0 through  the  relation 


(ACp) P = - 40X()  (17) 

Tlte  limits  of  integration  of  liquation  (16)  are  dependent  on  whether  the 
leading  edge  is  subsonic  or  supersonic.  Each  of  these  cases  will  be  considered 
separately  below. 
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Supersonic  Leading  Edge.  By  supersonic  lcac.1  inp  edge,  one  means  the  Mach 
number  normal  to  the  leading  edge  is  greater  than  one.  Referring  to  figure  4.  this 
implies  the  Mach  line  emanating  from  the  root  chord.  OA.  lies  behind  the  wing 
leading  edge.  If  the  Mach  line  intersects  the  wing  tip.  as  shown  in  figure  4.  there 
are  five  distinct  disturbance  regions  present  and  hence,  five  different  perturbation 
solutions.  It  the  Mach  line  intersects  the  wing  trailing  edge,  the  disturbance  created 
by  the  Mach  line  OA  impinging  on  the  wing  tip  is  eliminated  so  that  onl\  four 
perturbation  solutions  are  needed. 


REGION 


Figure  4.  Flat  Plate  Wing  Planform  With  Supersonic  Leading  and 
Trailing  Edges;  Mach  Line  Intersects  Wing  Tip 


The  generalized  formulas  for  the  lifting  pressure  distributions  on  a wing  in 
steady  roll  are  derived  in  Reference  21.  For  convenience,  the  final  equations  are 
repeated  here.  For  region  I.  the  (low  is  two-dimensional  and  the  resulting  lifting 
pressure  is 


(ACp)p  = 


(ACp)p 


4pm2  x0(m2o  - 1) 
i ~ p2  V^tm2  - 1)3/2 


(18) 


In  region  2.  the  total  lifting  pressure  is: 


<ACp)p  = 


(ACp)p 


4pm " x0  \ , , 1 + m2o 

= — z , rrr  \( I + lira)  cos  

2 7t)32V  (m2  - 1)3/-  | Lnv(l+o) 


-( 1 - nr  a)  cos 


i 


I - nr  a 
m(  l-o)  . 


(19) 


The  induced  pressure  caused  by  the  tip  Mach  line  l)F  is: 


( 

- 

- 

-niX| 

4 pm  l 

mxa  2 

b , 
-(nr  - 1) 
- 

0 + ya(2m  + 1 ) 

p,plj  ‘ 7rV„^(m2-l)3/2| 

-my  - 

L e 

cos 

mxa 

_ 


Jnv 


mya  \ 

(x  +0y  )(m  + 1 ). 


CO) 


I'he  total  lifting  pressure  in  region  3 is  then: 


(ACp)p  = ((AC'plpl,  + [(ACp)p  13 


(21) 


The  lifting  pressure  in  region  4 is  a combination  of  that  in  Regions  2 and  3.  Thus, 
in  Region  4: 


(ACp)p  = | ( AC'p )p  1 2 + [(ACp)p  1 3 


(22) 
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Again,  if  the  Mach  line  OA  intersects  the  wing  trailing  edge,  lapiations  18  through 
22  allow  one  to  determine  the  complete  lifting  pressure  distribution  over  the  wing 
surface.  If  the  Mach  line  OA  intersects  the  wing  tip,  another  perturbation  is  induced 
in  the  How  field.  The  total  pressure  differential  in  region  5 is  then: 


4pm  \ 

|,)p  = Tr/Wjm2-!)3'2  | 


nix,,  , b , 

+ in'  y + - (m“  + I ) 

P a 2 


cos  1 


mxa/j3  - ya(l  - 2m)  + b 
mx  /0  + ya  + b 


■ T 

-2m\  ya  ( m - 1 >(mxa//3  + iny;|  + b)  f 


(23) 


Subsonic  Leading  Edge,  If  the  Mach  number  normal  to  the  leading  edge  is  less 
than  one  (and  the  trailing  edge  Mach  number  is  greater  than  one),  there  are  only 
two  flow  regions  to  consider.  However,  the  solution  is  complicated  somewhat  by  the 
leading  edge  singularity  where  the  velocity  goes  like  l/x.22  Referring  to  Figure  5, 
the  pressure  differential  of  upper  and  lower  surfaces  in  Region  I is:22 


1 ODFG 

2 DFH 


Figure  5.  Flat  Plate  Wing  Planform  With  Subsonic  Leading  Edge  and 
Supersonic  Trailing  Edge 


13 


(24) 


2l<m)pnr  x0o 

(A(p)p  = 7T1T7 


Km)  = 


2(  1 - nr  ) 

(2  - nr  )H(m)  - nr  1( m) 


F(m)  = 


J 1 - k2  sin-  060 


F(m)  = j 


[ 1 - k’  sin’  0 


k = Jl  - nr 


Values  of  the  complete  elliptic  integrals  of  the  first  and  second  kind  (F(m)  and 
I (in),  respectively)  have  been  tabulated  and  appear  in  standard  mathematical 
handbooks. 

Region  2.  which  is  affected  by  the  wing  tip  Mach  cone,  has  a lifting  pressure 
differential  given  by: 

Spin  r r3mxo/0  + y(!  - 2m)- b/2(  I +m) 

<A('  )„  = — 1 — Jb/2-v  2 , , ■ ■■=  (25) 

it  L 3(l+m)J(l+m)(mx0/<3  + y0 ) . 


For  a How  with  a subsonic  leading  edge  and  supersonic  trailing  edge,  liquations  (24) 
and  (25)  determine  the  complete  pressure  distribution  due  to  a steady  roll.  p. 

The  local  roll  moment,  total  rolling  moment,  and  roll  damping  moment  for 
both  the  subsonic  and  supersonic  leading  edge  eases  can  then  be  determined  by 
liquations  (12).  (13)  anil  (14).  respectively. 


Transonic  Flow  (Mfh  < M,  < 1.2) 


There  are  currently  no  simple,  accurate  analytical  methods  available  for 

calculating  transonic  roll  damping.  With  the  increased  emphasis  on  transonic 
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aerodynamics  in  the  last  lew  years,  it  is  hoped  a simple  theoretical  method  will 
become  available  within  a reasonable  time  frame.  Until  that  time,  one  must  resort 
to  the  lull  Navier  Stokes  equations  of  motion  or  empirical  techniques.  Since  the 
former  approach  is  beyond  the  scope  ol  this  work,  the  empirical  approach  will  he 
followed. 


Without  a theoretical  model  to  calculate  transonic  roll  damping,  it  has  been 
current  practice  by  some  engineers* ' to  estimate  roll  damping  in  direct  proportion 
to  the  lift.  That  is 


(C 


" (C\,»'M  Mtb 


N„'m 


lb  (C 


N 'M  = \! 


fb 


(26) 


where  the  quantities  C\  > and  CN  are  total  configuration  values.  This  means  that  if 
the  roll  damping  is  known  at  say  = 0.8  and  ( N is  known  throughout  the 

transonic  speed  regime,  then  the  roll  damping  can  be  estimated  according  to 
I quation  (26).  This  is  the  procedure  used  in  the  present  analysis  for  calculating  roll 
damping. 

One  of  the  main  problems  inherent  in  calculating  transonic  aerodynamics  is  a 
need  to  account  for  thickness,  and  in  some  cases,  aei  ■ effects.  An  empirical 

technique  proposed  by  l.dmondson-^  to  correct  lor  thickness  and  aeroelastic  ettects 
on  rectangular  wings  is 


ft.  = 


1>  iheor  y 


(I  - t/c) 


,»2  AK  / 3 


(27) 


where  t/c  is  the  overall  thickness  to  chord  ratio  of  the  wing  planform.  However,  the 
present  values  of  K’v  , \i,cor>  *n  transonic  How  are  calculated  based  on  the  lilt 
curve  slope  which  has  already  accounted  tor  thickness.  In  addition,  intuitively,  « ne 
would  expect  a Mach  number  effect  to  be  present  in  1 quation  (27),  since  thickness 
is  most  important  in  transonic  flow,  for  these  two  reasons,  f quation  (2  ) was  not 
applied  in  the  present  analysis. 


I nterlereiuv  I ffects 

There  are  three  types  of  interference  effects  which  need  to  be  examined  in  the 
process  of  predicting  roll  damping.  These  are  the  body-fin.  tail  tin  to  tail  tin  and 
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canard  to  tail  interference  effects.  The  first  arises  due  to  the  presence  of  the  body, 
the  second  as  a result  of  tail  shed  vortices  impacting  another  tail  fin,  and  the  third 
as  a result  of  canard  shed  vortices  impacting  the  tail  tins. 

Using  slender  body  theory,  the  effect  of  the  body  has  been  calculated  by 
Adams  and  Dugan25  for  both  planar  and  cruciform  wing-body  combinations.  The 
results  of  their  calculations  are  presented  in  Figure  6 as  a function  of  the  parameter 


£ = rb/(rb+b/2) 

It  is  seen  that  there  is  little  difference  between  the  planar  and  cruciform  wing 
configuration  and  furthermore,  that  the  body  has  little  effect  on  the  roll  damping 
for  values  of  £ up  to  0.4.  Although  these  results  are  analytical,  they  have  been 
verified  qualitatively  by  the  experimental  results  presented  in  Reference  I . 
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Figure  6.  Effect  of  Body  Radius  on  Damping  in  Roll  and  Pitch  for 
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The  effect  of  the  number  of  fins  on  roll  damping  was  at  first  determined  using 
slender  body  theory.27  The  results  indicated  that  the  addition  of  fins  to  a missile 
adds  to  the  damping  in  roll  at  a decreasing  rate.  However,  slender  body  theory 
results  are  independent  of  Mach  number.  Moreover,  most  tactical  weapons  fiy  at 
Mach  numbers  such  that  the  nondimensjona!  spin  pfi  f/2Voo  is  so  low  that  tail  fin 
shed  vortices  do  not  affect  other  tail  fins.  Hence,  tail  fin  to  tail  fin  interference  for 
most  practical  cases  is  zero,  and  is  therefore  neglected  here. 

Canard-tail  interference  will  not  be  accounted  for  at  this  time.  Several  personnel 
are  currently  involved  in  estimating  induced  roll  and  it  is  hoped  that  this  form  of 
interference  can  be  incorporated  into  the  current  methodology  at  some  future  time. 


MAGNUS  KORIT:  AND  MOMENTS 

There  have  been  several  analytical  attempts  at  predicting  Magnus  (References  28 
through  31 1 characteristics  as  a function  of  the  main  variables  of  interest;  that  is. 
Mach  number.  Reynolds  number,  body  shape  and  various  properties  of  the  boundary 
layer  All  of  these  methods  have  given  some  success  for  given  conditions,  but.  as  of 
yet,  there  is  still  no  accurate  analytical  method  by  which  one  can  calculate  Magnus 
as  a function  of  Mach  number  and  body  shape  for  typical  projectile/missile 
ordnance.  The  latest  of  the  theoretical  attempts  to  predict  Magnus  is  due  to 

Vaughn.31  I'o  evaluate  the  capability  of  Vaughn's  method,  several  projectile 
configurations  were  considered  in  which  his  method  was  compared  with  experimental 
data  and  an  all-empirical  prediction  model  “Spinner”.4  figure  7 presents  the  results 
of  these  calculations  As  seen  in  the  figure,  the  empirical  model  is  better  for 

predicting  Magnus  than  the  theoretical  model  of  Vaughn.  Since  neither  approach 

does  very  well  in  general,  it  appears  there  is  still  much  analytical  work  to  be  done 

before  Magnus  forces  and  moments  on  projectile  and  missile  ordnance  can  be 
predicted. 


Reference  32  attempted  to  generalize  the  empirical  Spinner  approach  to 
incorporate  boattail  geometry  variations  in  the  prediction  scheme.  Graff3"  et  al 
found  that  one  could  correlate  the  magnus  force  and  moment  as  a linear  function 
of  the  Mangier  variable,  ??.  where  r?  is  defined  by; 


V 


5 r2dx 

i Mi  1 

r?cr  Vr  J 


(28) 
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which  for  conical  boattuils  integrates  to 


i?  = 


(29 ) 


An  empirical  scheme  was  thus  developed  which  included  boattail  shape  in  addition 
to  length  m the  methodology.  However,  when  the  method  was  compared  to  other 
prediction  techniques  for  existing  ordnance,  it  appeared  the  method  ot  Reference  22 
did  not  do  significantly  better  than  the  Spinner  technique.  Since  the  method  in 
Reference  32  was  only  for  angles  of  attack  up  to  1°,  the  Spinner  program  will  be 
used  to  estimate  magnus  effects  on  spin-stabilized  ordnance. 

lor  fin-stabilized  weapons  it  will  be  assumed  the  configuration  rolls  slow 
enough  so  magnus  effects  can  be  neglected. 


PITCH  DAMPING  MOMENT 

There  are  two  contributions  to  the  damping  in  pitch.  These  are  the  steady 
state  derivative  CM  (the  moment  coefficient  due  to  steady  pitching  velocity)  and 
the  unsteady  state  or  time  dependent  derivative  ( M (the  moment  coet t icient  due 
to  a constant  vertical  acceleration).  In  order  to  determine  the  total  damping  ot 
missile  configurations,  the  damping  of  individual  components  (i.e.,  body,  wings,  and 
canards)  is  calculated,  and  then  the  wing-body  interference  is  included. 

The  empirical  procedure  of  the  “Spinner”  program  is  employed  to  estimate  the 
body-alone  pitch  damping  parameter.  Thus,  the  quantity  CM()  + (m-  *or  ;l  particular 
body  is  related  to  the  damping  value  for  a dimensionally  similar  body  by  a function 
of  the  body  length,  boattail  length,  and  center  of  gravity. 

The  lilting  surfaces  of  missiles  (canard-body-tail  configurations)  and  fin-stabilized 
projectiles  (body-tail  configurations)  generally  represent  the  major  contribution  ot  the 
total  damping  in  pitch.  Therefore,  it  is  important  to  have  a good  estimate  ot  the 
wing-alone  pitch  damping. 

lor  subsonic  and  supersonic  flow  calculations,  once  again  a thin  wing  is 
assumed  so  that  small  perturbation  theory  is  applicable.  Also,  as  in  the  section  Roll 
Damping  Moment,  the  wing  is  assumed  to  have  no  camber  or  twist  and  negligible 
aeroelastic  effects.  In  transonic  How,  empirical  methods  are  applied  in  longitudinal 
damping  computations.  The  component  interference  etlects  are  then  approximated  by 
assuming  a slender  body  and  using  the  apparent-mass  method  ot  Bryson. 
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The  individual  procedures  used  for  calculating  the  wing  pitch  damping  moment 
coefficient  derivatives  in  all  flight  regimes  under  consideration  and  the  interference 
effects  are  presented  in  the  following  discussion. 

Subsonic  Flow  (M  < M,i 

°° fb 

The  stability  derivatives  C,  and  for  a thin  wing  in  subsonic  flow  are 

j 1 o M o 

determined  front  the  pressure  distribution  corresponding  to  a steady  pitching  rate  q. 
In  the  case  of  incompressible  How  the  pressure  variation  is  obtained  by  solving 
I -((iiation  (I),  the  small  perturbation  equation  for  three-dimensional  steady  flow,  with 
the  appropriate  boundary  condition.  For  steady  state  pitch  damping  calculations  the 
necessary  boundary  condition  is  found  by  setting  a(|  = p = 0 in  Equation  (5). 
Thus, 


*z0  = C1(X0  - Xref)/V 


( 30 ) 


Moreover,  the  pressure  loading  due  to  a constant  pitch  rate  is  realized  by  solving 
the  integral  equation  given  by  Equation  (II)  with  py /V^  replaced  with 

q( \0  - xr(,|)/Voo.  The  resulting  solution  is  the  wing  loading  ACp.  Compressible 
subsonic  calculations  are  based  upon  the  incompressible  solution  by  applying 
(iothert’s  extension  to  the  Prandtl-Glauert  rule,  which  is  given  along  with  the 
expression  for  the  pressure  coefficient  C„  in  Equations  (/)  and  (9),  respectively. 

I rom  the  chordwise  pressure  distribution,  the  sectional  moment  about  the  axis 
of  rotation  of  the  wing  due  to  q is  determined,  specifically. 


CM 
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i T1' 

1 x(  A(  „ ) dx 

c8  J H 4 

ret  x . - 


(31 ) 


where  the  quantity  (ACp)  represents  the  loading  due  to  a pitch  rate  q.  The  total 
pitching  moment  due  to  q on  the  entire  wing  is  given  by: 
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\ssuming  the  pitching  moment  is  a linear  (unction  o(  pitch  rate,  the  steady  state 
pitch  damping  moment  coefficient  derivative  is 


(33 » 


Many  of  the  recent  developments  in  the  field  of  unsteady  aerodynamics  (i.e.. 
References  34  and  35)  indicate  that  oscillatory  lifting  surface  theory  represents  the 
most  viable  approach  in  the  prediction  ol  unsteady  subsonic  stability  derivatives  tor 
arbitrary  wing  planforms.  However,  the  computational  effort  required  by  this  theory 
to  obtain  such  derivatives  precluded  the  possibility  of  being  able  to  readily 
incorporate  the  technique  into  the  present  subsonic  calculation  scheme  (steady  lifting 
surface  theory).  Two-dimensional  unsteady  methods  (i.e..  References  36  and  37)  tor 
estimating  . and  C’M  . , which  assume  small  amplitude  oscillations  and  correct  for 
finite  span  effects,  appear  to  be  inadequate  to  handle  basic  wing  planforms  (i.e., 
rectangular,  triangular,  tapered,  sweptback,  etc.)  for  a complete  aspect  ratio  range. 
Although  a simple  solution  for  C’M  . in  the  case  of  a triangular  wing  is  obtained  by 
Tobak  and  Lessing  in  Reference  38,  the  mathematical  intricacies  associated  with  the 
method,  which  reduces  an  unsteady  flow  problem  to  an  equivalent  steady  flow 
problem,  prohibit  similar  simple  solutions  for  other  wing  planforms. 


In  the  present  treatment  of  subsonic  How,  the  constant  vertical  acceleration 

derivative  (\,  . for  the  wine  alone  is  given  its  steady  flow  value  (CM.  = 0).  Since 
° ‘ <* 
the  primary  concern  here  is  the  prediction  of  dynamic  derivatives  for  tail-body  or 

wing-body-tail  configurations,  this  deficiency  does  not  seriously  impair  the 
preliminary  design  objective.  It  is  interesting  to  note  the  theoretical  predictions  of 
( M which  assume  CM  . = 0,  for  a delta  wing  (4R  = 2.4)  and  a sweptback 

and  tapered  wing  (A3  = 3)  that  are  presented  in  Figures  20  and  21.  These  figures 
are  discussed  in  the  results  section. 


If  only  small  amplitude  oscillations  are  considered  (i.e..  reduced  frequency 
cud  2V  < 1 ).  and  the  flow  over  the  bodv  is  assumed  to  be  in  phase  with  the 
flow  over  the  tail  (i.e.,  the  time  required  by  a change  in  body  vortex  strength  due 
to  a change  in  angle  of  attack  to  travel  to  the  tail  is  small  in  comparison  to  the 
period  of  oscillation),  it  is  probably  reasonable  to  assume  that  the  unsteady 
derivative  (\.  . is  approximately  zero  for  a tail-body  in  subsonic  flow.  However,  this 
is  not  the  situation  for  a canard-body-tail  configuration.  In  particular,  the  tail 
contribution  to  the  damping  parameter  CM  . resulting  from  the  downwash  lag 


phenomenon  should  he  considered.  Since  there  is  .1  finite  time  required  lor  the 
dovenwash  discharged  from  a forward  lifting  surface  to  reach  the  tail,  it  can  he 
shown'  that  the  tail  contribution  to  CM  is  given  by 
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where  the  subscript  T means  the  value  associated  with  the  tail,  and  xrcf  is  the 
distance  between  the  vehicle  center  of  gravity  and  the  center  of  pressure  of  the  tail. 
The  downwash  variation  with  respect  to  angle  of  attack  can  be  calculated  from  the 
tail  effectiveness  parameter.  That  is. 
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where  C is  the  lift  coefficient  for  the  negative  lift  at  the  tail  due  to  the 
downwash  induced  by  the  wing  shed  vortex,  (CL)BT  is  the  lift  coefficient  for  the 
body-tail  combination,  and  Ct  is  the  lift  coefficient  for  the  body  alone.  Details 
for  computing  these  quantities  are  given  in  Reference  7. 


Supersonic  Flow  > 1.2) 


The  pressure  loading  distributions  to  calculate  the  pitch  damping  CM  ( +CM. 
are  obtained  by  using  liquation  (17).  To  determine  the  perturbation  velocity 
potential  and  thus  0x  in  that  equation  the  unsteady  wave  equation,  which  is 
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atul  is  derived  from  linearized  small  disturbance  theory,  is  solved.  On  the  basis  of  a 
first-order  frequency  solution  to  the  above  equation,  an  expression  tor  the 
time-dependent  disturbance  velocity  potential  for  a wing  accelerating  downward  at  a 
positive  angle  of  attack  at  is  given  by'4 


</>  = 


fixing  the  same  analysis  it  can  also  be  shown  that  the  lilting  pressure  (evaluated  at 
time  t 0)  is 


,,  = 

where  0(  , and  < A('p ) ( , are  the  perturbation  velocity  potential  and  the  pressure 
loading,  respectively,  due  to  unit  (positive)  steady  pitching  velocity  about  the  y-axis 
and  where  <p  , and  ( ACp  )u  , are,  analogously,  the  perturbation  velocity  potential 

and  pressure  loading,  respectively,  due  to  unit  (positive)  angle  ol  attack.  I bus.  the 

unsteady  flow  problem  is  reduced  to  an  equivalent  steady  flow  problem,  and  it  now 
remains  to  determine  the  steady  quantities  just  mentioned  in  order  to  obtain  the 
desired  stability  derivatives  and  pitch  damping. 

1 rom  the  results  that  are  given  in  References  40  and  41,  the  potential  and 

corresponding  ACp  due  to  angle  of  attack  are  readily  determined  for  a thin  wing  in 
supersonic  flow.  In  particular,  Equation  ( I ).  which  is  the  steady  flow  version  ol 
I quation  (35).  is  solved  by  I ward’s  source  sheet  method4"  using  the  appropriate 
boundary  condition  ( <i> t = a).  Ihe  solution  is  given  in  I quation  (lb),  lor  a constant 
rate  of  pitch  the  solution  just  referred  to  is  again  applicable  with  the  associated 

boundary  condition  (0  = qx/V_  ).  As  it  has  been  indicated  previously,  the  limits  ol 

integration  of  I quation  (lb)  are  dependent  on  whether  (lie  leading  edge  of  a given 
wing  is  subsonic  or  supersonic. 
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Supersonic  Leading  Edge.  A description  of  the  supersonic  leading  edge  ease  is 
given  m the  Supersonic  Leading  Fdge  for  Roll  Damping  Moment.  The  generalized 

formulas  for  the  lifting  pressure  distributions  on  a wing  with  a constant  pitching 

velocity  or  a constant  vertical  acceleration  are  derived,  except  for  region  5.  in 

References  3‘>  and  40.  Referring  to  Figure  4.  the  How  in  region  I is 

two-dimensional,  and  the  resulting  lilting  pressure  is 
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In  region  2.  the  total  lifting  pressure  is 
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The  induced  pressure  caused  by  the  tip  Much  line  DF,  added  to  the  contribution  ol 
region  1,  yields  the  total  lifting  pressure  in  region  3: 
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file  total  pressure  differential  in  region  5 is 


S, . I, sonic  Leading  Edge.  For  a wing  with  * 
description  is  given  in  the  Subsonic  Leading 
Referring  to  Figure  5.  the  pressure  differentials 

region  1 are4  1 


subsonic  leading  edge  the  physical 
1 dge  for  Roll  Damping  Moment, 
of  the  upper  and  lower  surlaces  in 
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where 
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and  the  elliptic  integrals  I (ml  and  I (in)  are  defined  alter  I quation  (24). 


Region  2.  which  is  affected  by  the  wing  tip  Mach  cone,  lias  lifting  pressures 
given  by 
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l or  a llow  with  a subsonic  leading  edge  and  supersonic  trailing  edge,  liquations  (48) 
through  (51)  determine  the  complete  pressure  distribution  due  to  steady  pitching 
rate  q and  constant  vertical  acceleration  a (~z). 

The  local  pitching  moment  due  to  q and  the  total  pitching  moment  due  to  q 
are  calculated  using  liquations  (31)  through  (33).  The  same  quantities  due  to  a are 
obtained  using  liquations  (31)  through  (33)  by  replacing  q and  i)  with  a. 


Transonic  Flow  ( M)h  < M < 1.2) 

The  problem  of  no  simple,  accurate  analytical  method  to  estimate  dynamic 
aerodynamic  coefficient  derivatives  is  encountered  again.  As  in  previous  cases,  the 
only  viable  approach  consistent  with  the  present  computational  cost  effective 
objective  is  an  empirical  one. 

The  first  step  in  the  approximation  of  the  transonic  pitch  damping  is  to 
determine  the  value  at  the  force  break  Mach  number  by  applying  the 

aforementioned  subsonic  calculation  procedure.  Then,  employing  the  empirical 
method  discussed  in  Reference  7,  sweep,  Mach  number,  aspect  ratio,  and  thickness 
are  taken  into  account.  That  is 
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where  the  quantities  CM  + C’M  and  C’N  are  the  total  configuration  values. 

Moreover,  if  the  damping  in  pitch  is  known  at  the  nominal  force  break  Mach 
number  of  0.8,  and  CN  is  known  throughout  the  transonic  speed  regime,  an 

estimate  of  pitch  damping  in  transonic  flow  can  be  made  with  liquation  (52). 


Interference  bffects 

As  mentioned  earlier,  once  the  aerodynamic  loadings  on  the  components  of  a 
flight  configuration  due  to  a particular  kind  of  motion  are  estimated,  the 

interference  effects  are  then  determined  in  order  to  obtain  the  total  forces  and 

moments.  An  application  of  the  apparent-mass  method  of  Bryson’’  provides  an 

approximation  for  an  interference  damping  reduction  factor,  which  is  defined  as  the 
ratio  of  the  pitch  damping  of  a wing-body  combination  to  the  pitch  damping  erf  the 
total  wing  (i.e.,  the  exposed  wing  panels  extended  through  the  body).  The  variation 
of  this  factor  with  respect  to  the  nondimensionalized  body  radius  £,  which  is 

defined  in  the  Interference  I fleets  for  Roll  Damping  Moment,  is  presented  in 

I igure  tv  This  figure  shows  a decrease  in  the  damping  of  the  total  wing  as  a result 
of  the  presence  of  the  body.  Moreover,  the  damping  reduction  factor  reaches  a 
minimum  of  .75  at  £ equal  to  .71.  After  this  facto-  is  determined  for  a given 

configuration,  the  total  wing-body  derivatives  CL(j,  CM(),  CL.,  and  (\,  are  found 

by  adding  the  modified  wing  value  to  the  exposed  body  alone  contribution,  where 
the  exposed  body  is  the  section  of  the  body  beyond  the  apex  of  the  wing. 
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When  a configuration  includes  both  forward  and  aft  lifting  surface  panels,  the 
induced  damping  in  pitch  due  to  wing-tail  interference  requires  attention  also.  Tlte 
effects  of  the  wing  shed  vortices  on  the  tail  contribution  to  the  longitudinal 
damping  can  be  calculated  in  one  of  two  ways:1 

I.  If  there  is  direct  interaction  between  the  aft  lifting  surfaces  and  the 
vortices  shed  from  the  forward  lifting  surfaces  (bw/bT  < 1.5),  the  wing-tail 
interference  is  approximated  by 


where  (C.  ) ..  is  the  negative  lift  of  the  tail  due  to  wing-tail 

L a w ( v ) 

interference,  and  x , is  the  distance  between  the  vehicle  center  of  gravity 
and  the  center  of  pressure  of  the  tail.  The  sign  of  the  above  expression  is 
negative  if  there  is  pure  pitching  motion,  and  it  is  positive  if  there  is 
vertical  acceleration  motion. 

2.  If  there  is  no  such  direct  interaction  (b^  /bT  > 1.5),  the  wing-tail 

interference  is  neglected. 

Note,  only  one  trailing  vortex  is  assumed  for  each  forward  lifting  surface  panel. 

figure  8 summarizes  the  methods  used  for  computing  the  dynamic  derivatives. 
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Figure  8.  Methods  Used  to  Compute  Dynamic  Derivative 
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UKRINI  ST  ATI’S  AND  I I'TI  KI  II  I OK  I 


The  current  status  of  the  present  effort  to  develop  an  aerodynamic  prediction 
capability  for  tactical  weapons  is  shown  in  figure  9.  As  seen  in  (lie  figure, 
aerodynamics  can  be  estimated  up  to  low  supersonic  Mach  numbers  ( M „ ' 2.M  and 
angles  of  attack  to  20°  for  nonairbreathing  tactical  weapons. 
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Figure  9.  Current  Status  of  Aerodynamic  Prediction  for  Tactical  Weapons 


It  appears  that  future  tactical  weapons  will  be  designed  for  higher  speeds, 

incorporating  higher  angles  of  attack  for  improved  maneuverability.  To  obtain  the 

higher  speeds,  airbreathing  configurations  appear  to  be  one  ot  the  more  attractive 

methods.  The  main  thrust  of  the  future  effort  will  therefore  be  to  extend  the 

methodology  to  high  supersonic  Mach  numbers  for  angles  of  attack  to  180"  and  to 
also  consider  the  airbreathing  configurations. 


I wo  other  areas  are  also  of  interest  for  future  effort.  These  are  improving 
those  methodologies  which  have  been  outdated  in  the  past  several  years  by  new 
technology  and  optimization  of  the  computer  program  for  ease  of  use  when  the 
work  is  complete  It  is  envisioned  these  three  elforts  will  take  Irom  3 to  4 years, 
dep  tiding  on  the  funding  available  in  a given  year,  starting  in  fiscal  year  l‘> 
\ltcr  the  work  is  complete  and  documented,  the  computer  program  will  be  turned 
- i |.  ilu  lomputer  programming  division  at  NSWC  1)1  lor  maintenance,  updating. 
. i t i ii  suppling  to  other  parts  of  the  Dot),  along  with  their  contractors 
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STATU  AERODYNAMICS 

Three  cases  are  presented  to  show  comparison  of  the  present  static 
aerodynamics  methodology  with  experimental  data.  These  cases  consist  ol  a spin 
stabilized  projectile  (body-alone  configuration),  an  unguided  missile  (body-tail 
configuration),  and  a complicated  canard-body-tail  configuration.  These  cases  are 
sufficiently  different  and  complex  to  test  all  the  individual  theoretical  and  empirical 
procedures  of  I igures  2 and  3.  and  to  indicate  typical  results  to  be  expected  from 
the  prediction  of  static  forces  and  moments.  For  more  detailed  comparisons,  refer  to 
References  5 through  l>. 

Figure  10  gives  the  static  aerodynamic  coefficients  for  the  improved  5"/54 
projectile.  I'he  improved  round  has  a 2.75  caliber  nose  and  a 1.0  caliber  boattail 
with  a discarding  rotating  band.  Txcellent  agreement  with  experimental  data  is 
obtained  lor  the  drag  coefficient  throughout  the  entire  Mach  number  range.  Fair 
agreement  is  obtained  for  normal  force  coefficient  and  hence  pitching  moment  and 
center  of  pressure.  The  comparison  for  the  lifting  properties  is  Mach  number 
dependent:  in  the  low  supersonic  region  the  theory  is  consistently  about  10  percent 
low  on  normal  force,  whereas  at  intermediate  supersonic  speeds  it  compares  very 
well  with  experiment.  The  reason  for  this  behavior  is  the  failure  of  the  in  viscid 
theory  to  predict  afterbody  lift  correctly  at  low  supersonic  Mach  numbers.  At 
subsonic  and  transonic  Mach  numbers,  the  theory  does  about  as  well  as  could  be 
expected,  considering  the  amount  of  empirical  methodology  in  that  region. 

The  second  case  is  a 10-caliber  missile  with  clipped  delta  tail  tins.  The 
experimental  data  are  taken  from  Reference  43.  which  gives  the  static  aerodynamics 
for  OX  < M < 1.3.  Figure  I)  compares  theoretical  drag  coefficient,  normal  force 


Figure  10.  Comparison  of  Theory  and  Test  Data  for  Improved 

5754  Projectile 


coefficient  derivative,  and  center  of  pressure  with  experiment  as  a function  ot  Math 
number  and  for  a = I ’.  Recall  from  figure  3 that  for  > 1.2.  the  hit  and  mag 
(except  for  base  drag)  was  calculated  numerically,  whereas  for  0.8  < M < 1.2  the 
theory  consists  of  mostly  empirical  procedures.  I or  < 0.8.  the  wing  lilt  is 
calculated  analytically  but  most  other  force  components  are  computed  empirically, 
With  the  exception  of  the  normal  force  coefficient  slope  at  M_  = 0.8  and  0 85.  the 
theory  is  well  within  10  percent  of  experiment.  The  maximum  error  in  center  of 
pressure  for  this  configuration  is  5 percent  of  the  length  or  half  a caliber. 

The  final  example  chosen  is  a complex  canard-body-tail  configuration.  Fhe  body 
nose  is  60  percent  blunt  with  two  ogive  segments  and  a 0.7  caliber  boattail.  The 
canard  has  an  aspect  ratio  of  two  with  a sweepback  angle  of  15°.  Its  shape  consists 
of  a sharp  wedge  leading  edge  with  a constant  thickness  section  following.  The 
trailing  edge  is  truncated  parallel  to  the  leading  edge.  The  tail  has  an  aspect  ratio  of 
four  with  cylindrical  leading  and  trailing  edges.  The  tail  thickness  to  chord  ratio  also 
varies  along  the  span.  The  detailed  canard  and  wing  geometry  listed  above  is  not 
needed  in  calculating  lift,  but  it  must  be  known  for  drag  computations.  The  results 
of  the  calculations  for  this  configuration  are  shown  in  Figures  12A  and  I2B. 
Figure  12A  gives  the  normal  force  and  center  of  pressure  for  Mw  = 1.6  and  at 
various  angles  of  attack.  Four  curves  are  shown  in  the  figure:  canard-body-tail  with 

canards  deflected  up  by  10’,  canard-body-tail  with  no  canard  deflection,  body-tail, 

and  finally,  body  alone.  Several  points  are  worthy  of  note  in  this  figure.  First  of  all 
the  body  alone  solution  agrees  very  well  with  the  unpublished  experimental  data  up 
to  a = 16°.  Above  a = 16°,  the  theory  is  low  which  is  probably  due  to  not  taking 
into  account  Reynolds  number  effect  in  the  body  crossflow  drag  coefficient.  The 
next  point  is  that  for  this  configuration,  the  tail  lift  is  about  10  percent  too  high 
and  the  canard  lift  about  15  percent  too  low  so  that  the  total  lift  agrees  almost 
perfectly  with  the  experimental  data  up  to  the  point  where  stall  begins  to  occur 

(a  > 14').  This  in  turn  causes  the  center  of  pressure  to  be  up  to  approximately 

0.75  calibers  more  rearward  than  the  experimental  data  suggest.  It  is  suspected  that 
the  theory  being  high  for  the  high  aspect  ratio  tail  and  low  for  the  moderate  aspect 
ratio  canard  is  due  to  the  flowfield  interaction  effects  from  the  complex 
configuration  and  will  not,  in  general,  be  true  for  other  cases.  However,  it  does 
indicate  that  the  theory  can  be  used  quite  effectively  in  design,  even  for  quite 
complex  wing-body-tail  geometries.  The  final  point  to  be  emphasized  from 
I figure  I2A  is  the  fact  that  no  attempt  has  been  made  to  predict  stall 
characteristics.  As  seen  in  the  figure  lor  this  configuration,  stall  occurs  around 
a = 15°  at  M„  = 1.6.  However,  if  the  wing  thickness  or  freestream  Mai  h number 
is  changed,  the  stalling  angle  of  attack  will  also  change. 
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Figure  12A.  Normal  Forces  and  Center  of  Pressure  of  a Missile  Configuration; 

AH.  = 4,  /R  = 2,  M = 1.6 

t C ' ® 


34 


1.0  1.5 


Figure  12B.  Drag  of  a Missile  Configuration  and  Its  Components 


Hie  drag  characteristics  for  this  same  missile  are  shown  in  figure  I2B.  flu- 
drag  is  shown  as  a function  of  Mach  number  and  again  the  total  force  is  broken 
down  into  its  components:  body-alone,  body-tail,  and  canard-body-tail.  The 

body-alone  drag  is  acceptable  in  supersonic  and  subsonic  flow  but  is  unacceptable  in 
transonic  flow  where  the  empirical  nature  of  the  theory  does  not  account  for  nose 
bluntness  correctly.  Flic  wing-alone  drag  shown  at  the  bottom  figure,  includes  the 
increase  in  base  drag  due  to  tails.  This  causes  the  tail  drag  to  be  high  because  the 
theory  predicts  this  base  drag  increase  to  be  significantly  higher  than  the 
experimental  data  suggest.  However,  the  body-tail  drag  is  still  within  the  ± 10-percent 
category.  Finally,  the  canard  drag  shown  at  the  bottom  figure,  is  added  to  the 
body-tail  drag  and  the  overprediction  of  tail  drag  is  compensated  somewhat  by  the 
under  prediction  of  canard  drag. 


ROLL  DAMPING  MOMENT 
Comparison  With  Exact  Linear  Theory 

The  present  method  of  computing  roll  damping  moment  in  supersonic  How  is 
to  compute  the  pressure  coefficients  and  integrate  these  numerically  over  the  entire 
wing  planform  (assume  wing  extends  through  body).  This  is  as  opposed  to 
computing  Cv„  for  simple  planform  geometries  using  very  lengthy  closed  form 
solutions.  The  advantage  to  the  former  approach  is  that  if  one  is  interested  in 
pressure  coefficients  and  local  wing  loadings,  in  addition  to  roll  damping  coefficients, 
it  is  easier  to  numerically  integrate  the  known  pressure  coefficients  than  to  program 
the  lengthy  equations  for  CViJ. 

In  subsonic  How,  closed  form  solutions  are  not  possible  so  no  comparisons  with 
exact  theories  can  be  made.  However,  the  methodology  used  is  basically  that  of 
References  17  and  44,  with  minor  modifications,  and  extensive  checkout  has  been 
performed  in  those  references. 

Empirical  methods  are  used  for  the  damping  moment  calculations  in  transonic- 
now  so  that  comparison  with  experiment  will  have  to  suffice  for  check-out  of  the 
approach. 

Three  cases  are  considered  as  test  cases  to  compare  the  numerical  solutions 
with  closed  form  analytical  solutions  such  as  presented  in  References  21  and  22. 
Fhese  include  a wing  with  subsonic  leading  and  supersonic  trailing  edges  and  a wing 
with  supersonic  leading  and  trailing  edges  with  the  Mach  line  intersecting  the  tip 


2 b 


and  then  the  trailing  edge,  tacit  of  these  cases  is  sufficiently  different  so  as  to 

check  the  present  numerical  results  with  the  closed  form  analytical  solutions  these 
cases  are  presented  in  Figures  13,  14,  and  15.  The  results  are  shown  as  the  local 

normal  force  coefficient,  cN)(,  as  a function  of  the  position  along  the  winy 

semispan.  The  local  normal  force  coefficient  due  to  a rolling  rate  p is  defined  by 


„ >P 

Pb/2V„ 


As  expected,  the  numerical  results  duplicate  the  analytical  solutions  in  all  cases. 


ANALYTICAL 

NUMERICAL 


Figure  13.  Wing  With  Subsonic  Leading  and  Supersonic  Trailing  Edges 

(cN  versus  y/b/2) 


Comparison  Willi  Experiment 


Much  work  has  been  done  in  measuring  roll  damping,  both  with  rocket 
powered  techniques  in  the  early  1950’s  and  in  the  last  tew  years  with  the  wind 
tunnel.  Comparison  of  these  experimental  results  with  the  3-1)  thin  wing  theory  in 
supersonic  How  and  lifting  surface  theory  in  subsonic  flow  has,  in  general,  shown 
reasonable  agreement.  I lie  theory  typically  overestimates  the  actual  experimental 
results  on  wings  where  aeroelastic  effects  are  important,  particularly  in  the  transonic 
and  low  supersonic  speed  regimes. 

Three  examples  in  which  roll  damping  have  been  calculated  are  shown  in 
I igures  16,  IN.  and  19  Figure  16  is  a comparison  of  theory  and  experiment  for 
the  delta  wing  configuration  of  Reference  45.  As  seen  in  the  figure,  the  present 
methodology  gives  reasonable  agreement  for  wing  t/c  values  of  0.04.  Apparently, 
aeroelastic  effects  were  important  for  this  configuration  because  the  linear  theory 
was  higher  than  experiment  by  as  much  as  25  percent.  I'll  is  is  consistent  with  the 
results  reported  in  References  48  and  49  which  indicated  aeroelastic  effects  can 
reduce  roll  damping  by  up  to  25-30  percent,  particularly  in  transonic  flow. 

G EXPERIMENT  It  c 04  REF  45 
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Figure  16.  Comparison  of  Theory  and  Experiment  for  Roll  Damping  of  a 

Delta  Wing 
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Figure  17.  Body-Tail  Configurations 
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Figure  18.  Roll  Damping  Coefficient  for  Army  Navy  Finner 
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Figure  19.  Roll  Damping  of  Navy  Research  Model  Configuration 


A second  comparison  ot  experiment  and  the  present  methodology  is  shown  in 
f igure  IK  I his  ease  is  for  the  cruciform  wing-body  Army-Navy  Pinner,  which  is 
illustrated  in  1 igure  17.  therefore,  the  interference  theory  procedures  can  also  be 
employed.  Considerable  experimental  data  exists  for  this  case  at  supersonic  speeds. 
In  addition,  the  quaxi-two-dimensional  approaches  ol  Nicolaides48  and  Oberkampf49 
are  shown  lor  comparison.  The  theoretical  approaches  are  all  within  the  experimental 
accuracy,  lending  support  to  the  linearized  supersonic  flow  theory.  Since  there  were 
no  data  available  for  this  case  in  subsonic  flow,  the  present  theory  could  not  be 
compared  with  experiment. 

In  f igure  17  a research  development  model  of  a beamrider  projectile  is  shown 
in  addition  to  the  Army-Navy  Pinner.  Since  the  present  analysis  lor  supersonic  flow 
does  not  handle  nonstreamwise  tips,  the  tips  of  the  beamrider  are  assumed  to  be  in 
the  streamwise  direction.  Moreover,  the  assumed  tip  chords  puss  through  the 
midchord  points  of  the  actual  tip  chords,  and  thus  the  wing  surface  area  is 
maintained.  The  variation  in  roll  damping  with  Mach  number  is  given  for  this 
configuration  in  Figure  19.  There  are  no  available  damping  in  roll  data  for  the 
supersonic  region.  The  theoretical  prediction  in  this  region  compares  favorably  with 
an  estimate  based  upon  the  experimental  lift  and  center  of  pressure. 


Pi t ch  I ) amping  Moment 


Comparison  With  Exact  Linear  Theory.  For  the  reason  indicated  in  the  roll 
damping  section,  the  pitch  damping  in  supersonic  How  is  computed  by  integrating 
the  pressure  distributions  corresponding  to  the  two  types  of  motion  associated  with 
longitudinal  damping.  A comparison  is  made  in  Figures  20  through  25  between  the 
numerical  and  analytical  results  for  the  spanwise  variation  of  the  sectional  normal 
force  derivative  due  to  q motion  and  the  sectional  normal  force  derivative  due  to  a 
motion.  These  quantities  are  defined  in  the  following  way; 


where 


c is  the  mean  aerodynamic  chord. 
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Figure  20.  Wing  With  Subsonic  Leading  and  Supersonic  Trailing  Edges 

(cNq  versus  y/b/2) 
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Figure  21.  Wing  With  Subsonic  Leading  and  Supersonic  Trailing  Edges 

(cN  versus  y/b/2) 
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Figure  24.  Wing  With  Supersonic  Leading  and  Trailing  Edges;  Mach  Line 
Intercepts  Tip  (cNq  versus  y/b/2) 
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Figure  25.  Wing  With  Supersonic  Leading  and  Trailing  Edges;  Mach  Line 
Intercepts  Tip  ( cN  versus  y/b/2) 
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I'liori’  are  negligible  differences  between  the  analytical  and  numerical  results  in  the 
figures. 


The  comments  concerning  subsonic  flow  and  transonic  flow  that  are  made  in 
the  previous  section  are  also  applicable  in  the  case  of  pitch  damping. 


C omparison  Wit  h I- \ peri  men  t 

As  a consequence  of  the  reduction  in  longitudinal  damping  experienced  by  most 
flight  vehicles  traveling  at  transonic  speeds  and  undergoing  short  period  oscillations,  a 
number  of  experiments  were  conducted  m the  1 950’s  to  investigate  pitch  damping. 
Several  experimental  techniques  ti  e . wind  tunnel,  ground  launched  rocket  models, 
full  scale  flight  tests)  were  utilized  to  obtain  measured  values  of  the  damping  in 
pitch.  At  the  present  time  most  of  these  tests  are  performed  in  the  wind  tunnel 
using  either  the  free  or  forced  oscillation  technique.  Comparisons  made  between 
these  data  and  the  linear  prediction  methods  discussed  herein  have  shown  reasonable 
agreement.  That  is.  trends  are  generally  predicted  well,  and  magnitudes  are  usually 
predicted  within  roughly  ±20  percent. 

In  Figures  26 a and  26b  the  subsonic  variation  of  the  quantity  CM  +CM.  with 
Mach  number  for  two  delta  wings  which  have  aspect  ratios  of  two  and  four, 
respectively,  are  presented.  Although  there  is  substantial  scatter  in  the  data,  the 
steady  damping  term  CM  is  considered  to  provide  an  acceptable  estimate  of  the 
total  damping.  Also,  the  calculation  of  for  a sweptback  and  tapered  wing  of 

aspect  ratio  three  and  taper  ratio  0.6,  which  is  shown  in  Figure  27.  agrees  fairly 
well  with  experiment  in  the  subsonic  region.  However,  the  peaking  of  the  pitch 
damping  in  the  transonic  region  is  not  predicted  very  well.  This  is  essentially  due  to 
the  application  of  an  empirical  calculation  procedure  that  is  based  upon  linear 
theory  and  the  neglect  of  the  unsteady  term  C'M  . , which  becomes  more  important 
in  the  transonic  region. 

To  further  evaluate  the  present  methodology  for  determining  the  pitch  damping 
the  following  three  wing-body  combinations,  which  arc  shown  in  Figures  28  and  29. 
arc  considered:  (I)  delta  wing-body  < A3  = 2).  (2)  delta  wing-body  (A3  = 4). 
(3)  sweptback  and  tapered  wing-body  (A3  = 3).  Since  these  three  wing-body 
configurations  do  not  represent  actual  tactical  weapons  (i.e..  the  location  of  the  axis 
of  rotation  from  the  nose  of  the  body  in  terms  of  percent  body  length  is  too 
high),  the  body  contribution  to  the  total  pitch  damping  is  a theoretical  estimate 
rather  than  the  usual  empirical  one.  Test  data  for  CMo  + CM.  as  a function  of 
Mach  number  and  the  corresponding  theoretical  predictions  are  presented  for  these 
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Figure  26.  Pitch  Damping  of  Two  Delta  Wings  Oscillating  About  Axis 

Through  cR  / 2 
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DIMENSIONS  IN  CALIBERS 
I CALIBER  = 10  16cm 


Figure  29.  Delta  Wing-Body  and  Sweptback  and  Tapered  Wing-Body 


three  configurations  in  i igures  30  through  32.  for  all  three  cases  the  agreement 
between  theorv  and  experiment  is  fairly  good  in  the  supersonic  speed  regime.  I he 
subsonic  prediction  is  reasonable  in  all  cases  since  the  maximum  deviation  is  about 
22  percent.  In  case  ( I ) the  transonic  peak  is  underestimated,  and  the  predicted 
force  break  Mach  number  is  higher  than  that  depicted  by  the  data  (i.e..  theory  gives 
peak  at  M = 1.14  while  experiment  gives  peak  at  M = .08).  As  mentioned 
previously,  this  is  a result  of  the  empirical  nature  of  the  transonic  scheme  and  the 
assumption  that  (\,  is  approximately  zero.  Configurations  (2)  and  (3)  encounter  a 

a 

region  of  instability  in  the  transonic  region.  The  method  used  here  to  obtain  pitch 
damping  indicates  it  a configuration  is  stable  or  unstable,  which  is  sufficient  with 
respect  to  current  stated  objectives,  but  it  does  not  predict  the  magnitude  of  the 
instability. 

finally,  the  pitch  damping  of  the  Army-Navy  Finner  and  a research 
development  model  of  a beamrider  projectile  are  examined.  Since  there  are  no 
longitudinal  damping  data  available  through  the  total  flight  range  for  a given  center 
of  gravity  location  in  the  case  of  the  basic  Finner,  theory  and  experiment  are 
compared  for  two  different  center  of  gravity  locations,  i.e.,  x,.(i  = 6.1  calibers  and 
x(  ( = 6.92  calibers.  Good  agreement  is  exhibited  in  Figure  33a  for  the  supersonic 

regime.  In  Figure  33b  the  rapid  decrease  in  the  damping  at  transonic  speeds  is  not 
predicted.  In  the  case  of  the  beamrider  the  pitch  damping  is  given  in  Figure  34. 
There  is  a maximum  deviation  of  23  percent  between  experimental  data  and  theory. 
The  decrease  in  the  theoretical  + CM  . in  the  lower  supersonic  Mach  number 

range  is  due  to  the  extension  of  the  wing  through  the  body  in  the  interference 
computational  procedure.  That  is,  there  is  a reduction  in  wing  tip  effects  because 
the  intersection  of  the  tips  by  the  Mach  lines  emanating  from  the  apex  of  the  total 
wing  are  downstream  of  the  intersections  by  the  Mach  lines  originating  from  the 
apexes  of  the  exposed  wing  panels.  The  maximum  reduction  in  total  damping,  which 
occurs  at  M = 1.2,  is  approximately  6 percent. 

COMPUTER  PROGRAM 

A computer  program,  which  utilizes  the  theory  discussed  in  this  report  and 
References  (6)  and  (7),  has  been  developed  to  calculate  both  the  static  and  dynamic 
aerodynamics  for  tactical  weapons  to  Mach  number  3 and  angle  of  attack  of  15°. 
The  primary  inputs  required  by  this  program  to  determine  the  aerodynamics  of  a 
given  llight  configuration  are  geometrical  ones,  i.e.,  body  length,  variation  in  body 
shape,  wing  planform  shape  in  terms  of  physical  dimensions,  etc.  Important  features 
of  this  program  are  its  fast  execution  time  and  low  cost,  lor  example,  it  takes 
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Figure  30.  Pitch  Damping  of  Delta  Wing-Body  (/R 
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Figure  31.  Pitch  Damping  of  Delta  Wing-Body  (/R 
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MACH  NO  - Mqq 

Figure  34.  Pitch  Damping  of  Navy  Research  Model  Configuration 


about  5 min  on  the  CDC  b70()  machine  at  a cost  of  about  $50  to  determine  the 
aerodynamic  coefficients  at  10  Mach  numbers  for  a configuration  such  as  the 
beamrider.  A detailed  description  of  the  computer  program,  along  with  several 
example  cases  and  a fortran  listing  of  the  program,  are  given  in  Part  II  of  the 
present  report. 


CONCLUSIONS  AND  RECOMMENDATIONS 


1.  A general  method  has  been  developed  consisting  of  several  theoretical  and 
empirical  procedures  to  calculate  lift,  drag,  pitching  moment,  magnus  moment,  roll 
damping  moment  and  pitch  damping  moment  on  wing-body-tail  configurations  from 
Mach  number  zero  to  three  and  angles  of  attack  up  to  stall. 

2.  Comparison  of  the  overall  methodology  with  experiment  for  several  configurations 
indicates  that  accuracies  of  ±10  percent  can  be  obtained  for  static  force  coefficients 
of  most  configurations.  Dynamic  force  and  moment  coefficients  are  within 
25  percent  of  experimental  data  for  most  cases.  This  is  at  a cost  of  less  than  $10 
per  Mach  number  for  the  entire  set  of  aerodynamic  coefficients. 

3.  If  roll  damping  is  desired  within  ±10  percent  for  most  configurations,  thickness 
and  aeroelastic  effects  should  be  accounted  for. 

4.  There  is  still  no  good  method  to  compute  magnus  forces  and  moments  on 
typical  ordnance. 

5.  The  prediction  methodology  in  transonic  How  is  almost  entirely  empirical  due  to 
the  current  state-of-the-art  in  simple  (but  yet  accurate)  analytical  methods  for  use  in 
flow  computations.  It  is  thus  recommended  that  work  be  continued  in  this  Mach 
number  regime  with  reasonable  emphasis  on  simple  techniques. 

6.  Due  to  increasing  performance  requirements  for  tactical  weapons,  it  now  appears 
appropriate  to  extend  the  present  design  methodology  to  incorporate  higher  Mach 
number,  higher  angle  of  attack,  and  airbreathing  configurations. 
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Co 


CM(1  + (M;. 


c'm. 


^ N 


CN 
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C 
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nh 


CN 
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Aspect  ratio. 

Wing  span  (does  not  include  body  radius) 

Chord  length  at  any  point  along  span. 

Mean  aerodynamic  chord 
Zero  lilt  drag  coefficient 

Rolling  moment  coefficient  of  wing  plan  form. 

Local  rolling  moment  coefficient. 

Roll  damping  moment  coefficient  |Cg/(p8rc(-/2V- )|. 

Pitching  moment  coefficient  measured  about  nose  tip  (positive 
up). 

Pitching  moment  coefficient  derivative  (d('M/da). 

Total  damping  in  pitch. 

Steady  state  damping  in  pitch  |(CM)  / (qCrcf/2V<jo ). 

Unsteady  damping  in  pitch  |(('M  )^/(atrcf/2Voo  )| . 

Normal  force  coefficient. 

Normal  force  coefficient  derivative  (dCN/dcr). 

Magnus  moment  coefficient  derivative. 

Magnus  force  coefficient  derivative. 

Sectional  normal  force  coefficient  derivative  due  to  p motion. 
Sectional  normal  force  coefficient  derivative  due  to  ()  motion. 


nose 


cN  Sectional  normal  force  coefficient  derivative  due  to  a motion. 

Pressure  coefficient . 

tA(\,l  Difference  in  pressure  coefficients  of  upper  and  lower  surfaces  of  wing 

planform. 

(A(p)p  Difference  in  pressure  coefficients  of  upper  and  lower  surfaces  of  wing 

planform  due  to  a rolling  velocity,  p. 

<ACpl()  Difference  in  pressure  coefficients  of  upper  and  lower  surfaces  of  wing 

planform  due  to  a pitch  rate  q. 

(M'pf  Difference  in  pressure  coefficients  of  upper  and  lower  surfaces  of  wing 

planform  due  to  a constant  vertical  acceleration  a. 

d Body  diameter. 

8 f Reference  length  (body  diameter  for  wing-body  configuration). 

liBT  Boattail  length. 

M Mach  number, 

m (3cot  A , . 

n Number  of  tail  fins, 

p Roll  rate, 

q Pitch  rate, 

q'  Dynamic  pressure. 

rb  Body  Radius. 

rn  Body  nose  radius. 

S Wing  planform  area. 

u.v.w  Perturbation  velocities  in  xQ,  y(),  z()  directions,  respectively. 
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franst'ormed  coordinates. 


Rectangular  coordinate  system  with  x at  nose  tip.  y out  right  wing, 
and  / positive  up.  If  configuration  consists  of  wing  only,  x begins  at 
wing  root  chord. 

Coordinates  of  wing  tip. 

Angle  of  attack. 

^ M2  - 1 . 

Canard  deflection  (degrees). 

Howarth-Mangler  variable. 

Boattail  angle. 

Sweepback  angle  of  a wing  generator  (i  = 1, 2,3.4)  with  i = I.  the 
wing  leading  edge,  and  i = 4 the  wing  trailing  edge. 

Taper  ratio. 

tan  A,  (y/x). 


Frequency  of  oscillation. 
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